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ABSTRACT 

The renormalization of general theories with inter-family mixing of Dirac and/or Majorana 
fermions is studied at the one-loop electroweak order. The phenomenological significance 
of the mixing-matrix renormalization is discussed, within the context of models based on 
the SU(2)i®U(l)y gauge group. The effect of radiative neutrino masses present in these 
models is naturally taken into account in this formulation. As an example, charged-lepton 
universality in pion decays is investigated in the heavy-neutrino limit. Non-decoupling 
heavy-neutrino effects induced by mixing renormalization are found to considerably affect 
the predictions in these new-physics scenarios. 
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1 Introduction 



Mixing effects have played a crucial role in the understanding of various aspects 
of particle phenomenology. The most well-known examples established by experiment 
are the K°-K° and B°-B° mixings [|IJ, which originate from the mixing of quarks ||. 
Furthermore, several astrophysical problems, including the solar-neutrino-deficit puzzle, 
may be explained by assuming that there is also mixing in the lepton sector leading to 
neutrino oscillations ||. In the Standard Model (SM), the mixing between the photon and 
the Z boson is implemented at the loop level, and precision test are becoming sensitive to 
this effect. In addition, many of the proposed extensions of the SM predict the possibility 
of mixings between scalar, fermion, and vector fields. In all these cases, the origin of the 
mixing is related to the rotation between weak and mass eigenstates. 

In order that a quantum theory yields precise quantitative predictions, it must be 
renormalizable. The renormalization of the SM has been established j|] and elaborated 

a long time ago. It was noticed [0 that the Cabibbo-Kobayashi-Maskawa (CKM) || 
matrix must be included in the renormalization programme as well. However, this effect 
has been found to be insignificant in the SM f9j. The reason is that the mass differences 
between the down-type quarks are small compared to the electroweak scale, so that the 
CKM matrix can effectively be taken to be diagonal. The situation should be very different 
in new-physics scenarios with large inter-family mixings. An attractive solution to the 
problem of the smallness in mass of the known neutrinos can arise from certain SO(10) 
grand unified models |ITJ] and/or E 6 superstring- inspired theories |TTJ, which predict see- 



saw-type neutrino mass matrices with large Dirac components fL2 |. 

The renormalization of mixing effects in extensions of the SM has not yet been ad- 
dressed in the literature. In this paper, we shall take the first step in this direction by 



elaborating the renormalization of general theories with Dirac and/or Majorana neu- 



trinos. Our formulation will naturally include radiative neutrino-mass contributions 
These considerations will affect a number of low-energy and LEP1/SLC electroweak ob- 
servables which have been utilized to establish bounds on the parameter space of these 
models. For illustration, we shall estimate the size of charged-lepton non-universality in 
pion decays. Specifically, we shall consider the heavy-neutrino limit of the observable 
R w = T(ir + -> e + u)/Y(Ti + -> fi+v). 

This paper is organized as follows. In Section 2, we shall introduce the formalism 
of mass, wave-function, and mixing-matrix renormalization in general models with Dirac 
and/or Majorana neutrinos and derive the corresponding counterterm (CT) Lagrangian. 
In Section 3, we shall determine the CT's for the Dirac case in the on-shell renormalization 
scheme. Special attention is paid to the balance of the numbers of CT's and renormal- 
ization conditions. These considerations are extended to the case of Majorana neutrinos 
in Section 4. In Section 5, we shall discuss the nature of the mixing matrices and their 
relationships in the framework of SU(2)^®U(l)y theories. The renormalization of these 
mixing matrices will then be performed in Section 6. We shall see that those relationships 
among the mixing matrices that are enforced by the unitarity of the theory carry over to 
the one-loop level, while additional identities related to other symmetries are violated by 
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ultraviolet divergences. As an application, in Section 7, we shall quantitatively analyze 
R n = r(7r + — > e + z/)/r(7r + — > y^v) in the heavy-neutrino limit. Our conclusions will be 
summarized in Section 8. 



2 One-loop renormalization 



We consider a general model with Nf fermions, / = (f\, . . . , f'N f ), which may be of 
Dirac and/or Majorana type. As usual, we define the left- and right-handed components as 
fi,R = Pl,rI, where Pl,r = (l=F7 5 )/2 are the corresponding chirality projection operators. 
We denote the weak eigenstates and their mass matrix by a prime. Bare parameters carry 
the superscript "0" . The bare kinetic Lagrangian in the weak basis is given by 

£L = M2 + M ~ fS^fn - fSM*«f2, (2.1) 

where 1 is the unity matrix in flavour space and M'° is a complex, non-diagonal mass 
matrix. M'° can always be diagonalized with the help of two unitary matrices, U® and U R , 
via 

M° = U L M'°U°J , (2.2) 

where M° is a non-negative, diagonal mass matrix. Consequently, the mass eigenstates are 
given by 

fl,R = Ul R f L \. (2.3) 

Next, we study the renormalization of this theory. For this end, we write the un- 
renormalized parameters and fields in terms of renormalized ones and CT's as 

U%R = U L , R + 6U LtR , (2.4) 
M'° = M' + SM', (2.5) 
M° = M + 5M, (2.6) 

f' L °R = z2gf' LtR = (l + l -5Z'^ fa, (2.7) 
fin = Zl%f LiR = (l + Uz L ^ f L , R . (2.8) 

The renormalization is arranged so that the basic structure of the theory is preserved 
through the order considered. U^r are still unitary, so that 

Ul,rSUI r + 6U LiR Ul iR = (2.9) 

is satisfied up to terms of 0(5Ul R ). Similarly to M°, M is diagonal with non-negative 
eigenvalues, and Eqs. (|2.2|) and ( |2.3| ) are mapped into 

M = U L M'U R , 
h,R = U LtR f' LjR , (2.10) 
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1/2 

respectively. It is sufficient to consider 5Ul,r, oM, and Z LR , since the renormalization 
constants connected with the weak eigenstates are determined by 

5M' = UlM8U R + Ul8MU R + 5UlMU R , (2.11) 
<2 = UI\ R Z]! 2 R U L ^ (2.12) 

Equation (|2.12|) is nicely illustrated in the following diagrams: 

r 7 ot r 7 ot 
fro l f o r/o y R no 

JL i JL JR i JR. 

Z'l l2 \ ]Z 1 ' 2 Z lX /\ ]ZT • (2-13) 

f , Ul f t i Ur r 

JL * JL JR > JR 

In the following, we shall work in the mass basis. In order to find the appropriate CT 
Lagrangian, we eliminate the bare masses and fields in Eq. (|2.1|), 



£«„ = ihZ^Z]! 2 fif L + if R Z^Z R /2 tf R - f L Zl /2 \M + 8M)Z R /2 f R 

-f R Z R /2 \M + 5M)Z l L /2 f L . (2.14) 

Up to higher-order terms, we then have 

f-'kin = {-'kin ~\~ $£>kim (2-15) 

where L^ n is the renormalized Lagrangian and 

&c*» = \h{sz L + sz^)9SL + \h{sz R + sz^)9f R (2.16) 

- f L (SM + l -M5Z R + hsZ L ^M)f R - f R (5M + l -M5Z L + ±6Z**M)f L 

is the CT Lagrangian. 

The mixing-matrix renormalization is only important if unitary matrices different 
from unity appear in the couplings, such as the CKM matrix in the SM charged- current 
interaction or the mixing matrices which give rise to flavour-changing neutral currents in 
new-physics scenarios. To elucidate this point, we consider the charged- current interaction 
in a SU(2)x,®U(l)y model with right-handed neutrinos. In this model, the bare Lagrangian 
for the interaction of the W boson with the charged leptons, l iy and the neutrinos, rz i? is 
given by 

JC° W = -l=(W;fj^ + h.c, (2.17) 

where g is the SU(2)i coupling constant and 

Jl = J£17m»l = ^°tX- (2-18) 
Here, B° is a CKM-type matrix in the lepton sector defined as 

B° = U° L ' l U° L ' n \ (2.19) 
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where and U^ n are the bare unitary matrices that participate in the diagonalization of 
the charged-lepton and neutrino mass matrices, respectively. Notice that the renormaliza- 
tion of the weak coupling g and the H^-boson wave-function are universal and determined 
by other renormalization conditions ||, so that we may separately consider the renormal- 
ization of the charged current J M . Substituting Eq. ( |2.4j ) in Eq. ( |2.19| ), we find 



B + 5B, 



(2.20) 



where 



5B = BUl5Ul ] + 5U l L U^B . 



From Eq. ( |2.9| ) we know that U^SUf and 8U l L U£ are antihermitean. Using Eqs. 
( p.20p ) we can write Eq. Q2.18| ) as J° = J M + <5J M , where 



2.21) 
81) and 



<5J M = l Llfl ^B(5Z n > L + 5Z n ^) + ^(5Z l > l + 5Z w )b 



+B 
+ 



5Z n,Lf i + U2SU L 



~ A (5Z^ 
-{6Z l ' L -5Z l ^)+5U l L U 1 ^ 



B\n L . 



(2.22) 



Here, we have decomposed 5Z n,L and 5Z l,L into hermitean and antihermitean parts and 
collected all antihermitean CT's within square brackets. Obviously, the presence of SB 
is indispensable to cancel the ultraviolet (UV) divergences contained in the antihermitean 
parts of SZ n ' L and SZ l,L . This will also be illustrated in Section 7, where we shall calculate 
the one-loop correction induced in by massive neutrinos. We may fix SU^ and 5U l L by 
requiring that the square brackets in Eq. ( j2.22j ) vanish. This leads to 



L -8Z nM )U n L , 
L -5Z IM )U 1 L . 



(2.23) 



A similar condition was proposed in Ref. in connection with the CKM mixing of the 
SM. In order to enforce the UV finiteness of physical observables, it would be sufficient 
to replace the parentheses in Eq. ( |2.23| ) with their UV-divergent parts, evaluated at some 
renormalization scale fi. This would correspond to the MS renormalization prescription. 
By the same token, we might add finite, antihermitean matrices, c™- and q -, to the terms 
contained within the parentheses of Eq. (|2.23| ) and fix them by imposing some additional 
renormalization conditions. This would not affect the hermitean parts of SZ n,L and SZ l,L , 
and in particular their diagonal elements, which are arranged so that the fermion propaga- 
tors have unit residues. Although the use of Eq. ( |2.23| ) is not compelling, this prescription 
seems natural and we shall adopt it in the remainder of this paper. A detailed study of the 
implications of general mixing-matrix renormalization schemes will be given elsewhere. 
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3 Dirac case 



In the following, we shall study the mass and wave-function renormalizations in a 
general theory involving the mixing of Nf Dirac fermions. We denote the Dirac fermions 
by fi, with i = 1, . . . , Nf. There are two sources of imaginary contributions to the bare 
amplitudes. They can arise either from the possibility of on-shell cuts through the loop 
amplitudes (absorptive parts) or from complex mixing parameters (CKM-type couplings). 
Because of the hermiticity of the bare and renormalized Lagrangians, the CT Lagrangian 
must be hermitean, too. Therefore, only the dispersive parts of the one-loop amplitudes 
can participate in the renormalization procedure. Consequently, we shall only consider 
the dispersive parts of the two-point functions in the following. Furthermore, it will be 
understood that complex conjugation acts on the coupling constants in the flavour space. 

We start by considering the unrenormalized fj — > f] transition amplitude. Its most 
general form in compliance with hermiticity is[] 

=4Pl^M 2 )+ + P ^g(9 2 ) + PR^*(g 2 ), (3.1) 



where 

The renormalized counterparts will be denoted by a hat. As per construction, they will 



^(? 2 ) = ^V), Sg(g 2 ) = S J t(g 2 ). (3.2) 



satisfy relations similar to Eqs. (|3.1| ) and (|3.2| 



From the CT Lagrangian ( |2.16| ), we read off the relations between the bare and 



renormalized transition amplitudes, viz. 

^M 2 ) = Eg(g 2 ) + i((5ZS + ^ 
m<f) = T$tf) + l(6Z* + 6Z> 
Sg(g 2 ) = Sg(g 2 ) - l -{ mi 5Z^+m 3 5Zf:) - S^Sm,. (3.3) 

Next, we evaluate the renormalization constants by imposing the on-shell renormalization 
conditions on the renormalized transition amplitudes. Specifically, the non-diagonal ele- 
ments of SZfj and 5Z^ are determined by requiring that S^-(^) be diagonal if the external 
lines are put on their mass shells, while the diagonal elements are fixed in such a way that 
the residues of the renormalized propagators are equal to unity, i.e., 

tijiflu^q) = 0, (3.4) 
ik(q)EM = 0, (3.5) 

' -%{<l)Ui{q) = 0, (3.6) 



rrii 



iii(g)£«(0- = 0. (3.7) 

4 -mi 



*This form generalizes the one used in Ref. M, which is specific for the SM. 
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Conditions (§2§)-(g77|) imply that 



m^(m 2 ) + ±?*{m 2 ) 


= o, 


(3.8) 


^Sg(mJ) + Sg(mf) 


= o, 


(3.9) 


miSgK) + SgK 2 ) 


= o, 


(3.10) 


m,SgK) + E£*(m 2 ) 


= o, 


(3.11) 


+ s2K a ) + arn^K?) + sf (m 2 )) 








= o, 


(3.12) 



where S'(g 2 ) 



8Z, 



dX(g 2 ) / 'dq 2 . For % ^ j, we obtain from Eqs. 
2 



u)p:',jiinj) + 77iiTWjE^(mjf) + ///,>;'](///]) + ///,>:],"(///] 



-(EH 



8Z, 



R 



1 1 

mf — wta 
1 j 



mirrijE^im 2 ) + m 2 £f*(m 2 ) + m,Sf(m?) + m&%* (m*) 



3 V 1 



3 "tj ' 



3* 



(3.13) 



In the diagonal case % = j, the number of renormalization constants to be determined may 
be reduced by exploiting the following symmetry present in the CT Lagrangian Q2.14j ): 



7-1/2 



J9i 7 l/2 



„tfli 7 V2 



(3.14) 



where are real phases. In this way, we may, e.g., arrange for all 5Z~ to be real. Employing 
Eqs. (|3.8j )- (|3.12| ), we then obtain 



5m; 



1 r 



mi 



ZE(m 2 )-ZZ*(m?. 



—rrii 



Rli 



— m, 



1 

—rrij 
2 1 



*#K) + K) 
1 
2 



£f (m 2 ) + Ef'K 

Sf(m 2 )4 - 



mi 



(3.15) 
(3.16) 

(3.17) 



At this stage, it is instructive to investigate the mass-degenerate limit of two fermions, 
fi and fj. We observe that 5Z^ and bZf- are singular in the limit m ; — > m^. However, due 
to condition ( p.23| ), only the hermitean parts of SZfj and bZf- occur in the renormalization 
of physical observables. For mj = rrij, we find 



->:£(m 2 )-m 2 



Sj'K) + £g'(mj) 



•m, 



Sg'(m 2 ) + Eg*'(m 2 ; 



(3.18) 



m 2 



Sg'H) + £g'(mj) 



•m, 



Sg'(m 2 ) + Eg*'(m 2 ; 



(3.19) 
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which are indeed finite. For i = j, Eqs. ( |3.18| ) and ( |3.19| ) coincide with the hermitean parts 
of Eqs. ( |3.15|) and (|3.16| ), respectively, as they should. 

Finally, we verify that the number of renormalization conditions equals the number of 
CT's. To be specific, in a model with Nf fermions, there are ANf real conditions for all i — j 
(Eqs. (f3.8|) ( pTTT] ) collapse to three real conditions), and ANf(Nf — 1) real conditions for 
all i ^ j (namely Eqs. (|3.^ )- (|3~Tl| )). Thus, the total number of renormalization conditions 
is ANf. On the other hand, counting the number of independent real CT's, we have Nf 
mass CT's and ANj wave-function renormalization constants (namely 5Z^, SZ^, and their 
complex conjugates). From the number of CT's one has to subtract the Nf phases, 8i, 
in Eq. ( |3.14| ), which can be used, e.g., to render 8Z^ real. As a result, the total number 
of independent real CT's is ANj ||, which is equal to the one of the real renormalization 
conditions. 



4 Majorana case 

In this section, we shall study the renormalization of a general model with Nf Majo- 
rana neutrinos. In a way, this is a generalization of the Dirac case considered in the previous 
section, since a Dirac fermion may always be represented as a pair of mass-degenerate Ma- 
jorana neutrinos. By the same token, it is possible to describe the mixing of Dirac and 
Majorana neutrinos. The kinetic Lagrangian of the model is given by 

£L = \ (ifl ml + HI fifn ~ f°LM°fa ~ /l^VE) • (4-1) 
The bare and renormalized fermion fields satisfy the Majorana constraint, 

fl = (f° R f, h = (f R f, (4.2) 

where C stands for the charge-conjugation operation. Using Eq. ( |2.8| ), we thus find that 

Zl /2 = Z^ 2 *. (4.3) 

Now, we can express the bare Lagrangian (|4.1| ) in terms of renormalized quantities as 

if L Z l l 2] Z)! 2 0f L + i f R Z 1 L /2T Z 1 L /2 * 0f R - f L Z]! 2 \M + 6M)Zl /2 *f R 
-f R Z 1 L /2T (M + 5M)Z 1 L /2 f L \. (4.4) 

Decomposing this as C° kin = Ckin + 8£kin, we obtain the CT Lagrangian, 



S£ kin = -f L (5Z L + 5Z L ^ 0f L + -f R (5Z L * + 5Z LT ) 0f R (4.5) 
~Jl{SM + ^MSZ L * + l -5Z L ^M)f R - l -f R (5M + l -M5Z L + ^5Z LT M)f L . 



S 



The most general form of the fj — *> fi transition amplitude between fermionic Majo- 
rana states reads 



(4.6) 



where we have made use of the relations 



^M 2 ) = £jT(g 2 ), W 



(4.7) 



which follow from the Majorana condition wL 



From the CT Lagrangian (|4.5| ), we may read off the relationships between Sf (g 2 ^ 



T.fUq 2 ), and their renormalized counterparts, viz. 



1 / 



(4.8) 
(4.9) 



Imposing the on-shell renormalization conditions ( |3.4| 
constants for i ^ j, 



we obtain the renormalization 



mf — m, 
1 j 



and those for z = j, 



+ 



-E^(m 2 ) - 2m 2 S^'(m 2 
1 



III; 



iM/, 



™ 2 ) + Ejf *'« 



2mi 



i*f(m 2 ) 



E^(m 2 ) + I 



(4.10) 



(4.11) 
(4.12) 



A special situation arises if a Majorana neutrino is massless at tree level. In contrast to 
the SM Dirac case, this does not necessarily imply that the mass CT Srrii in Eq. ( |4.12 ) 
vanishes. In general, it will be positive and finite, i.e., the Majorana neutrino receives a 
mass via loop effects. Various mechanisms for generating radiative neutrino masses have 



been suggested in the literature [njTSfl; they are naturally implemented in our formulation 

Again, we ma? 
finite in the limit m. 



Again, we may verify that the physically relevant hermitean part of 5Z^ is indeed 



rrij, 



-4K)- 

+£fV 2 ; 



m ■ 



4V 2 ) + srK 



rrij 



Moreover, for i = j, we recover the hermitean part of Eq. (|4.11|) , 



-S^m 2 ) 



2m 2 E^'(m 2 ) 



ITLi 



K 2 ) + Ef*' k 2 ; 



(4.13) 



(4.14) 
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These observations reassure us of the self-consistency of our formalism. 

Let us finally count the number of independent renormalization conditions and CT's 
in our iV/-Majorana-neutrino model. For i ^ j, we have only two independent complex 
equations or four real conditions resulting from Eqs. (|3.8|) - (|3.11|) , due to the Majorana 
constraint (|4.2|) . This gives 2Nf(Nf — 1), where we only consider the cases % > j, so 
as to avoid double-counting due to the hermiticity of the renormalization conditions. In 
addition, there are 3Nf real relations for the diagonal transitions. So, we obtain 2Nj + 
Nf independent real renormalization conditions in total. On the other hand, taking the 
Majorana constraint ( |4.3|) into account, we count the same number of independent CT's. 
In fact, there are 2Nj wave-function renormalizations, SZ^ and 8Z^*, and Nf mass CT's. 



5 Mixing matrices in SU(2)^(8)U(l)y theories 

A minimal, renormalizable extension of the SM that can naturally accommodate 
heavy Majorana neutrinos is a model based on the SU(2)x,®U(l)y gauge group, in which 
lepton-number-violating AL = 2 operators have been introduced in the Yukawa sector by 
the inclusion of a number of Nr isosinglet neutrinos. The latter are sometimes called right- 
handed neutrinos because they are blind under SU(2)l. Here, we adopt the conventions 
of Ref. |T5[] and denote the bare isosinglet weak eigenstates by i/£ (with i = 1, . . . , Nr). 
Furthermore, we assume a number of Nq weak isodoublets. The quark sector of this model 
is similar to that of the minimal SM. The bare Yukawa Lagrangian that describes the 
neutrino sector reads 



C 



/ = + h-c, (5.1) 
where M' 0,u is a complex, symmetric mass matrix of the from 

m "- = ( ™V < ) • (5 - 2) 

It can always be diagonalized through the unitary transformation 

U°>» T M'°< v U°< v = M '" . (5.3) 

The nonnegative, diagonal matrix M°> u contains the bare neutrino-mass eigenvalues. The 
corresponding mass eigenstates are given by 

/ ,/0 \ N G +N R / ,Q C \ N G +N R 



,/OC 



E U°f*n%, .Jo = E U^n%. (5.4) 



Here, the first Nc mass eigenstates, z/j = Ui (i — 1, . . . , Nq), are identified with the ordinary 
light neutrinos (if N G = 3), and the remaining N R states, Ni = n i+Na (i = 1, . . . ,N R ), 
are the new neutral leptons predicted by the model. These neutral leptons are the heavy 



10 



Majorana neutrinos, which should be heavier than the Z boson, as they have escaped 
detection in production experiments at LEP1/SLC. The diagonalization of the charged- 
lepton mass matrix proceeds as outlined in Section 2. 

In our minimal model, quantum mixing effects enter via the interactions of the Ma- 
jorana neutrinos, rii, and charged leptons, k, with the intermediate bosons, W ± and Z, 
and the Higgs boson, H . In the 't Hooft-Feynman gauge, the bare Lagrangians of these 
interactions are 1TH 



N G N G +N R 



W 



7. 



V2 l=1 j=l 

n N G N G +N R 

° (G-r e e + h. c . 



4c?„ 



(5.5) 
(5.6) 
(5.7) 



T^r E 175 K + mJJSfeCS + i(mj - m?)9mCg. In? , (5.8) 



Afe+^fl 



N G +N R 



#° E ( m ° + ™°)&eCg. + z 75 (m° - mJ)3mC; n° , (5.9) 



where G ± and G° are the charged and neutral Higgs-Kibble ghosts, respectively, g is the 
SU(2)l coupling constant, c w is the cosine of the weak mixing angle, and B and C are 
Nc x (Nq + Nr) and (Nq + Nr) x (Nq + Nr) mixing matrices, respectively. The bare 
matrices are defined as 



4 



a 



N G 

E 

fc=i 

E u° k fu, 



E^W 
=i 

iV G 



kj 



(5.10) 
(5.11) 



fc=i 



Note that C^- is hermitean. Comparing Eqs. fl2.19|) and ( 5.10 ), we may identify 



yO,l 



U 



0,1 
L > 



U' 



U 



0,nT 



(5.12) 



Furthermore, B° and C° satisfy a number of identities, which will turn out to be crucial to 
warranty the renormalizability of our model, namely |T5|rE ,r7 



N G +N R 

E B hk B hk = $hhi 
k=l 



(5.13) 
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N G +Nr 

E = (5.14) 

k=l 

N G +N R 

E BlCl l = Bl (5.15) 
fe=i 

N G 

EA°X— Cg- (5-16) 
In addition, there are identities involving the Major ana-neutrino masses, 

N G +N R N G +N R N G +N R 

E <C° k C% = 0, E n»M = 0, E <KA = ^ (5.17) 

fe=l fc=l k=l 

These relations signify the presence of lepton- number- violating interactions, e.g., in the 



possible neutrinoless double-beta decay of a nucleus [18 



6 Renormalization of the mixing matrices 

Inserting Eq. (|2.23|) into Eq. Q2.21 ), we obtain a closed expression for the CT matrix 



of B, 

5B = -(5Z l ' L - 5Z l ' L ^)B --B(5Z n ' L - 5Z nM ). (6.1 



4V j 4 

By analogy, writing C& = Cy + SCij and requiring that Eq. Q5.16Q also holds true for the 
renormalized quantities, i.e., that 

N G 

Y.BfiBu = (6.2) 
i=i 

is correct to one loop, we find 

6C = ^6Z n ' L - 5Z nM ) C-hj (5Z n ' L - 5Z nM ). (6.3) 

With the help of Eqs. (|6.1|) and (|6.3| ), we may verify that the renormalized versions of 
Eqs. (|5TTJ)-(|^TJ), 



N G +N R N G +N R N G +N R 

E BhkB* 2 k = ^htej E CikCkj = Cij, E Bi k Cki = B u , (6.4) 
fe=i fc=i fc=i 

are valid at the one-loop level. On the other hand, the relations in Eq. ( |5.17| ) become UV 
divergent if we replace the bare parameters with their renormalized counterparts. This may 
be attributed to the fact that, in contrast to Eqs. (|5.13| )- (|5.16|) , the relations in Eq. (|5.17|) 
are not enforced by unitarity. 
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7 The observable R 



In order to illustrate the phenomenological significance of mixing-matrix renormaliza- 
tion, we shall now study possible violations of charged-current universality in the leptonic 
pion decays tt + — > l + v within the Majorana-neutrino mixing models described in Section 5. 
To simplify matters, we shall consider the limit where the novel Majorana neutrinos are 
much heavier than the intermediate bosons. We may then exploit the Goldstone-boson 
equivalence theorem |l9j formulated in the 't Hooft-Feynman gauge and include in the 
loop amplitudes only the massive Higgs boson, H, the massless Goldstone bosons, G^ and 
G°, and the heavy Majorana neutrinos, iVj. In this way, we may extract the leading elec- 
troweak radiative corrections of 0{Gpm 2 Ni ) and 0{GpMfj), which occur because Majorana 
neutrinos and the Higgs boson do not decouple in the high- mass limit. 

The decay 7r + (P) — > l + (k)u(p) is described by the parton-model process ud — > l + v in 



connection with the hadronic matrix element (0|c/(0)7 m 75m(0)|7t + (P)) = /V-P^, where f w is 
the pion decay constant. The corresponding vector-current matrix element is taken to be 
zero. The tree-level transition-matrix element reads 



mi 



M w 



mi 



-u„(p)P R vi(k), 



(7.1) 



where a is the fine-structure constant, V is the CKM matrix, and cr„ 



M 2 W 



/M 



In a full electroweak one-loop analysis, one needs to include the corrections related to the 
W-boson propagator, the Wud and Wlu vertices, and the udlu box as well as the renor- 
malizations of M^y, a, s w , the external quark and lepton fields, and the mixing matrices, 
V and B. For the time being, we ignore photon bremsstrahlung and corrections due to 
strong interactions. Neglecting the pion, quark, and external-lepton masses, the corrected 
matrix element takes the form 



r = T (l + 8™ v + O + % + %, 



(7.2) 



where % and T b are the one-particle-irreducible vertex and box amplitudes, respectively, 
and we distinguish between universal and lepton-flavour-dependent CT contributions. To 
simplify the notation, we include the propagator corrections in the universal CT. Specifi- 
cally, we have 



zuniv 
Oct 



n 



WW 



(0) - SM W 



M w 



2^ 

e 



6 s,. 



v.. 



ud 



SV ud + y( SZ d*d V udi + V u^ Z u lU 



SBl + ~-{6Z L n : v B\ ni + BIJZ^ 



(7.3) 
(7.4) 



where summation over the quark and lepton flavours u$, di, k, and rij is implied and [Q 



6M W 



neU ww (M^) , 
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5e 
e 

Ss w 



-n' AA (o) - 



s w n ZA (o) 



M 2 



1 c 2 

2 si 



n 



WW 



(M, 



z 

2 ■» 



nzz(Ml) 

Mi 



(7.5) 



The transverse gauge-boson vacuum-polarization contributions induced by the leptons of 
our SU(2)x,(8)U(l)y model may be found in Ref. |I7|]. The calculation considerably simplifies 
in the heavy-neutrino limit, 3> My/. According to the equivalence theorem, it is then 
sufficient to compute the vertex and lepton- mixing diagrams depicted in Fig. 1, while 
7ft = 0. Furthermore, we may project out the relevant vertex form factor by 



(7.6) 



The present knowledge of the radiative corrections to r(7r 
nicely summarized by Marciano and Sirlin 
uncertainties cancel in the ratio R n = r(7r 



20 1 



» l + v) in the SM has been 
The short-distance corrections and most 



-z/)/r(vr- 



Apart from helicity- 



suppressed terms of 0[(a/7r) (m^/m 2 ,) ln(m 2 /m 2 )], where the typical hadronic mass scale 
m p is used as a demarcation between short- and long-distance loop corrections, only QED 



corrections survive. The latter consist of a pointlike-pion contribution and a structure- 
dependent bremsstrahlung correction, which is suppressed by m\jm^ p . The leading term 
of the pointlike-pion contribution is given by —(3a/2Tr)ln(m 2 l /m 2 ,) ~ —3.7% and may be 
summed via the renormalization group; it makes up the bulk of the SM correction to R n . 
According to Ref. 



the current SM prediction is 



R, 



SM 



1.2352 ±0.0005) x 10" 



(7.7) 



where the error is mainly due to model dependence. As we shall see in the following, the 
presence of heavy Majorana neutrinos might be manifested by a significant shift in the the- 
oretical prediction of R n relative to R^ M . In turn, confrontation of the modified prediction 
with experiment, together with similar analyses for other low-energy and LEP1/SLC ob- 
servables, will allow one to improve the constraints on the parameter space of the Majorana- 
neutrino models under consideration. In our approximation, the one-loop-corrected value 
of Rn in the Majorana- neutrino scenario is given by 



R K = R° 7r {l + C), 



(7.f 



where R® emerges from the SM tree- level result by multiplication with \B eu \ 2 ^B^ 2 and 



(7.9) 



The universal contribution S^ w 



has dropped out. Analytic results for 5 l v and 5 l ct are listed 
in the Appendix. We note that S^ w and 8 V + 5 ct are separately UV finite. The significance 
of mixing-matrix renormalization becomes apparent by observing that 5 l v + 5 l ct contains 
the UV-divergent contribution from 5B* V . We may refine Eq. ( [7.8|) by including the SM 

20(1. For this end, we write R® = (\B ei 



radiative corrections of Ref 
r sm is given by Eq _ fly^ 



)R% M , where 
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For clarity, we shall henceforth consider a very simple realization of a SU(2)i®U(l)y 
model with Majorana neutrinos, which is of phenomenological interest and displays the 
essential features of mixing renormalization. We shall assume that the SM is extended by 
two heavy Majorana neutrinos, N\ and N2, which only mix with the leptons of one family, 
the first one, say. The other two families are assumed to be standard. In the notation of 
Section 5, this corresponds to Ng = 1 and Nr = 2. The relevant mixing-matrix elements 
are determined by solving Eqs. (|5.13[) -( p.l7|) , with the result that 



\B P 



1 



B, 



eNi 



B, 



,1/2 



eN 2 



1 P 



-1/4 R 



eNi, 



(7.10) 



C 



P 



1 + pV2 



c 



N 2 N 2 



-1/2 C 



N-i N-i 



c 



N X N 2 



-c 



N 2 Ni 



-1/4 C 



N1N1 



where p 



rn 



N 2 



jvnr N and s v £ measures the degree of light-heavy neutrino mixing |22 |. 



Without loss of generality, we may assume that p > 1. 



As mentioned above, we set 
This pattern of mixing may be 
motivated by the non-observation of the decay p —>■ e'j or the absence of p-e conversion in 
nuclei. Furthermore, a global analysis of low-energy data gives the following upper limits 



0. Then, Eq. Q becomes C = 2$te{5 e v + 6, 



ct I 



< 0.010 



< 0.0020 



at the 95% confidence level. We should note that electroweak corrections, mainly those 
originating from the SM, are taken into account only for a limited number of low-energy ob- 
servables, such as the muon decay width, in the evaluation of the bounds given in Eq. (|7. 1 1|) . 
Therefore, the so-derived bounds should be considered to be of tree-level accuracy. 

The results for S l v and 5 l ct given in the Appendix are valid for Mh arbitrary as long 
as m^jm^ 3> Mw- If we are interested in the non-decoupling effects of 0{Gftti 2 n ), we 
may put Mh = 0. In this limit, we have 



C 



G F {s v £ 



m 



N 2 



4vr 2 V2(l + Jpf 



f( P ) + (1 - y/fl)f 



(7.12) 



where 



/(*) 



(7.13) 



If 1 \ 11 
- 1 - — In 1 -x\ . 

2 V x 2 ) 1 1 4 2x 

Notice that, to the order considered, we may introduce Gf in Eq. ( ffljl) by substituting 
Gf = {j^oil \f2s 2 w M^ v )\ the adjustment proportional to Ar will only appear at the two- 
loop order because the tree-level expression for does not contain a. The correction C 
is throughout negative, takes on the simple form 



C 



3G F (s" L - 



m 



at, 



for p 



1Qti 2 V2 

1, and exhibits the limiting behaviour 



C 



G F {s v £Ym 



Ni 



1 - 



1 



VP, 



In p 



11 

y 



6^/p 



O 



(7.14) 



(7.15) 
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for p ^> 1. It is interesting to observe that C is of 0(Gj?m^ ), i.e., it scales with the 
squared mass of the lighter neutrino. 

We are now in a position to estimate how the value of s v f extracted from the 
measurement of R n is affected by the inclusion of loop effects due to heavy Majorana 
neutrinos. If we neglect these loop effects, we obtain the value (s£ e ) 2 by equating 
j^exp _ jjSMq_ (s£ e ) 2 ]. On the other hand, including these effects, we extract the corrected 
value {s u £f = (s u L e ) 2 + A(s£) 2 from R e * p = R s n M [l - (s£) 2 ](l + C), where C is evaluated 
with s v £. Consequently, the relative shift in (s£ e ) 2 is 

^-'-W'-C fe-lV (7.16) 



In the degenerate case, = mjv 2 = mjv, we then have 



(sl) 2 V 1TeV 
Inserting in Eq. ( [7.171 ) the perturbative upper limit |24 



16%xfeY. (7.17) 



s>w~m ^lTeV, (7.18) 

we find that {s v £) 2 is reduced by 16%. 

With the help of Eq. ( |A.2| ), we shall now investigate the dependence of 
-A(s v L e ) 2 /(s u L e ) 2 w (1 - R K /R^ M )/(s u L e ) 2 on m^, m^, and M H . We first consider the 
degenerate case, with p — 1. In Fig. 2, we show — A(s£ e ) 2 /(s^ e ) 2 as a function of m^r for 
M# = 100, 300, 500, and 1000 GeV, assuming (a) (s v L e f = 0.01 and (b) 0.007. We observe 
that, at low values of m^, the relative reduction of (s^ ) 2 strongly depends on M H ; it in- 
creases by almost one order of magnitude as Mu runs from 100 GeV to 1 TeV. Comparing 
Figs. 2(a) and (b), we see that, for in the TeV range, the relative reduction of (s 1 ^) 2 
is slightly smaller for the lower value of (s^ e ) 2 . The high-m^r regime is well described by 
Eq. (|7.14|) . Next, we study the non-degenerate case for Mh = 200 GeV. In Fig. 3, we dis- 
play the p dependence of — A(s£ e ) 2 /(s£ e ) 2 for m Nl = 0.5, 1, 2, and 5 GeV, again assuming 
(a) (s£) 2 = 0.01 and (b) 0.007. We see that, in compliance with Eq. (|7TT5|) , -A(s£ e ) 2 /(s£) 2 
grows quadratically with for p fixed, while it grows logarithmically with p for 
fixed. From Fig. 3(a), we read off a 20% effect for mjVi = 5 TeV and = 50 TeV. On 
the experimental side, the measurements at the Tri-University Meson Facility (TRIUMF) 
p5| and the Paul Scherrer Institute (PSI) [^6| yielded 

^(TRIUMF) = (1.2265 ± 0.0034(stat) ± 0.0044(syst)) x 10~ 4 , 

R e n xp (PSI) = (1.2346 ± 0.0035(stat) ± 0.0036(syst)) x 10~ 4 , (7.19) 

respectively. This represents a remarkable reduction in error, by a factor of 3, relative to 
the previous value ^7|, R exv = (1.218 ± 0.014) x 10 -4 . It is therefore reasonable to expect 



that, in the near future, experiments will become sensitive to 10% effects in (sj 
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8 Conclusions 



The renormalization of general theories with inter-family mixing between fermionic 
Dirac and/or Majorana states was studied to one loop in the electroweak on-shell scheme. 
Special attention was paid to the renormalization of the mixing matrix |7j, which plays a 
central role in such theories. Similarly to the renormalization prescription for the CKM 
matrix of the SM proposed in Ref. ||, we adjusted the mixing-matrix CT's in such a way 
that they precisely cancel the antihermitean parts of the wave-function renormalization 
constants. Our formulation naturally takes possible radiative-neutrino-mass contributions 
into account. 

The phenomenological implications of mixing renormalization for Majorana-neutrino 



mass models with large SU(2)i-breaking Dirac components [16| were analyzed in Section 7. 
In such scenarios, low-energy observables may receive sizeable corrections due to the non- 
decoupling of heavy neutrinos. As an example, the electroweak corrections to the observable 
R n were estimated in the heavy-neutrino limit. It was found that they may reduce the tree- 
level values of R n predicted in these models by up to 0.2%. This has to be contrasted with 



the present theoretical error in the SM prediction, which is ±0.04% |20|| . Future experiments 



to be performed at TRIUMF and PSI may be sensitive to such new-physics phenomena. 

Finally, we wish to emphasize that our formalism for the mixing renormalization of 
fermionic states may straightforwardly be extended to theories which involve the mixing of 
scalar or vector particles, such as the mixing of scalar quarks in supersymmetric theories. 
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A Analytic expressions 



In this appendix, we list the lepton transition amplitudes as well as the corrections 5 l v 
and 5 l ct defined in Eqs. (|7.4j ) and (|7.6| ), respectively, to one loop in the Majorana- neutrino 
models of Section 5 assuming mjsr i 3> Myy and Mg arbitrary. We adopt the Passarino- 
Veltman [^] conventions for the standard one-loop integrals in dimensional regularization, 
implemented with the Minkowskian metric, g^ u = diag(l, —1, —1, —1), as in Appendix A 
of Ref. ||. 

The charged- lepton and Majorana- neutrino mixing amplitudes depicted in Figs. 1(b) 
and (c) are found to be 



a 



8ns 2 w M 2 

+(m k P L + m h P R )m 2 k B (p 2 , m 2 k , M\ 



^{m 2 k P L + m h m h P fl ) B l (p 2 } m 2 k} M ] 



2 ) 
w) 



2 s 



a 



tfP L (miC* k + m k Cu?j (m k C kj + mjC* kj 



+ tfP R (miC ik + m k C* k ) (m k C* kj + m 3 C kj )] [SiQo 2 , m 2 , M|) + B^p 2 , m 2 k , M 2 H 
+ P L m k (miC ik + m k C*,?j (m k C kj + m,-Cjy) + P R m k (miC* k + m k C ik ) 

x (m fc C^. + m.Ai)] [5 (p 2 , m 2 , M|) - 5 (p 2 , m 2 , )] }, (A.l) 



respectively. Here and in the following, the Majorana indices rij are abbreviated by i, and it 
is summed over the indices of the heavy Majorana neutrinos. For the Wlv vertex correction 
due to Fig. 1(a) and the lepton-flavour-dependent CT contribution to T(7r + — > l + v) we find 



'ct 



a 



SttsIM 2 Cum * 



a 



C 24 (0, 0, 0, M|, m 2 , 0) + C 24 (0, 0, 0, 0, m 2 , 0) 



327TS2M 2 , 

+ |C;/m 2 



Cum 2 



B 1 (0,mlM 2 H ) + 3B 1 (0,mlQ) 



BiK 2 , m 2 , M 2 H ) + B 1 (ml m 2 , 0) - ^(0, m% M 2 H ) - B^O, m 2 , 0) 



-B (ml m 2 , M 2 H ) + B (m 2 , m 2 , 0) + B (0, m 2 , M 2 H ) - B (0 } m 2 , 0) 



+Cf j m mj [B 1 (mlm 2 ,Ml) + B 1 (mlm 2 j ,0) - -^(B (m 2 ,m 2 , M 2 H ) 



-B (ml m 2 , 0) - 50(0, m 2 , M 2 H ) + B o (0, m 2 , 0) 



(A.2) 
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Figure Captions 

Fig. 1: Feynman graphs contributing to the observable = T(n + — > e + v)/T{it + — > 
/i + z/) in the heavy-neutrino limit: (a) Wlv vertex corrections, (b) charged- 
lepton transition amplitudes, (c) neutrino transition amplitudes. 

Fig. 2: -A(s£) 7(s£ e ) 2 « (1 - i^/i?°)/(s£) 2 as a function of m N = m Nl = m N2 for 
M H = 100, 300, 500, 1000 GeV, assuming (s v £f = and (a) (s£ e ) 2 = 0.01 
or (b) (s£) 2 = 0.007. 

Fig. 3: -A(s£)7(s£) 2 « (1 - R n /R° 7r )/(s u L e ) 2 as a function of p = m 2 N Jm 2 Nl for 
mjVl = 0.5, 1, 2, 5 TeV, assuming M# = 200 GeV, (s^) 2 = 0, and (a) 
{s v £f = 0.01 or (b) (s£) 2 = 0.007. 
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